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, We present first-principles band calculations as well as structural optimization of the orthorhom- 

^\ ' bic KiCeo polymerized phase. We found three-dimensional inter-fuUerene bonding/anti-bonding 

0^ \ characters consisting of ti„ molecular orbitals with p-like symmetry in the conduction bands. The 

. formation of the four-membered ring connecting the fuUerenes lifts the px band from the continuous 

^ ' Py-Pz band, where the z-axis is parallel to the chain. The asymmetry between px and py may play 

U , an important role in binding the chains with the proposed rotational configuration ffni. 

Q. 

^ • 31.15.Ar, 61.46. +w, 61.48. -l-c, 61.50.-f, 71.20.-b 

^ ■ 

The recent discoveries of new alkali-doped fuUerene 

I ■ polymers § || § Q H show diversity in the ffeome- 

•7^ tries of the interconnections of fuUerenes. An impressive 
^y-j ' property of this polymer family is that the charge state of 
' ' the fuUerene seems to govern the polymer geometry. The 
AiCeo (A=K, Rb) g] and the NasRbCgo polymers 
Ch ' lH Hi have one-dimensional interconnections, while the 
""I , Na4C6o polymer has two-dimensional interconnections. 
Even in the case of the one-dimensional polymers, we can 
find a remarkable difference between the AiCeo (A=K, 
^ • Rb) and the Na2Rb'C6o polymers; the former is known to 
' have two covalent bonds between the fuUerenes , while 
^ ■ the latter is deduced from its larger inter-molecular sepa- 
O ' ration Q to have one covalent bond. These dramatic 
differences in bonding which are controlled by the charge 
states give rise to the question as how the excess electrons 
T— I on the fuUerenes affect the polymer geometry. 
^ ■ In this paper, we present a first-principles electronic 
I structure study of the ortho-KCeo polymer to elucidate 
how the electronic states of the conduction bands are re- 
, lated to its one-dimensional polymer geometry. Although 
. electronic band structure calculations of AiCeo have been 
— . ' reported by two other groups [f), to our knowledge, 
^\ this is the first quantitatively reliable calculation on the 
system with a realistic atomic configuration. The calcu- 
^ , lation method is based on the density functional theory 
^ • within the local density approximation (DFT-LDA) Q. 



A plane-wave basis expansion and optimized pseudopo- 
tentials in separable form were adopted in this simula- 
tion. 



I 

O - 

Q . Structural optimization is performed on the KCeo 
^ ' polymer, which has an orthorhombic unit cell containing 
. ^ , two Ceo molecules and two K atoms. The cell param- 
' eters are fixed at the experimental values [||. Initially, 
5^ ' the Cgo molecules are placed at the body-centered posi- 
tions, in which the atomic configuration of an isolated Ih 
Cgo molecule is used. We choose the initial orientational 
configuration of the fuUerenes as the one proposed by 
Stephens et al. ||l[], i.e., the unit cell is taken to be sim- 
ple orthorhombic. The structural optimization is then 
performed without any symmetry constraint. To obtain 
as precise a geometry as possible, we have checked the 
convergency of the geometry as a function of the cut- 



1 



off energy for plane waves ranging from 30Ry to 80Ry 
with uniform sampling of either one (F) or four A:-points. 
The geometry is relaxed until all the forces acting on 
the atoms become less than 5 x 10"'* [Hartree/au]. Co- 
valent bonds spontaneously appeared between the Ceo 
molecules during the optimization process, resulting in a 
one-dimensional chain structure. 

The calculated bond lengths are displayed in TABLE 1 
along with the experimental values. The calculated val- 
ues are consistent with the experimental values within 
the limit of error, with the largest discrepancy occurring 
in the Cl-Cl intra-fuUerene distance. Refinement of the 
experimental value (1.74 A) is desired, since it seems quite 
large compared with the typical C-C bond lengths known 
thus far. 

Let us discuss the geometrical character of the KCeo 
polymer from the viewpoint of the coupling between the 
conduction electrons and the deformation of the Ceo 
molecule. It would be natural to speculate that the char- 
acter of the LUMO state of isolated Cgo would be re- 
flected in the conduction band. Therefore, we shall note 
the symmetry of the LUMO state of isolated Ih Ceo be- 
fore explaining the resultant structure. 

The LUMO of the It Ceo molecule are threefold- 
degenerate tiu states with p-like symmetry. The coor- 
dinate axes for these p-like {tiy) states can be mapped 
on the three orthogonal C2 rotational pivots on the Ih 
Ceo molecule. 

If we now go back to the calculated polymer structure, 
we note that the chain direction is parallel to one of the 
C2 rotational pivots. Therefore, if we take the chain di- 
rection as the z-axis of the txu states, the other two co- 
ordinate axes are automatically determined: one (x-axis) 
is on the plane which includes the four-membered ring of 
the inter-molecular connection and the other (y-axis) is 
normal to the above two. We then measure the diame- 
ters of the polymerized fullerene along these coordinate 
axes. Interestingly, even the two diameters normal to the 
chain are significantly different from each other: they are 
6.98A and 6.83A along the x- and y-axes, respectively. As 
is expected, the diameter along the z-axis is elongated to 
7.48A. 

We emphasize here that this ellipsoidal distortion may 
be plausibly associated with the Jahn- Teller effect; how- 
ever, as will be clarified later in this paper, the relative 
positions of the ii„ levels in the polymer are instead de- 
termined by the crucial change in the local electronic 
states around the four-membered ring connecting the 
fullerenes as well as the strength of the inter-molecular 
interaction. 

The electronic states of this orthorhombic polymer 
crystal have been investigated using the first-principles 
band calculation method. The householder diagonaliza- 
tion method with a plane- wave expansion is used to main- 
tain the symmetry of the eigenstates as precisely as pos- 
sible. The cutoff energy for the plane waves is reduced 
to 25Ry to save on the computational cost. The num- 
ber of basis functions is then reduced to 25000. We note 
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here that the convergence of the band structure has been 
confirmed for both diamond and potassium crystals with 
this cutoff energy, though it is inadequate for the struc- 
tural optimization. Because the diagonalization of a ma- 
trix with a dimension of 25000 requires a huge amount 
of CPU time even for the newer parallel supercomputers, 
we did not perform fully a self-consistent diagonalization 
calculation but instead obtained the eigenvalues and the 
eigenstates through the following two steps. First, the 
effective potential for the Kohn-Sham equation was cal- 
culated self-consistently by means of conjugate gradient 
minimization of the total energy, in which 75 uniform 
A:-points in the first Brillouin zone were sampled. We 
then diagonalized the matrix calculated from the given 
effective potential. 

The calculated band structure is shown in FIG. 2. Al- 
though the band dispersion along the F-H line (parallel to 
the chain) is slightly larger than that along the F-M line 
(normal to the chain), the band structure indicates essen- 
tially three-dimensional character. This is qualitatively 
consistent with previous calculations |^ , although the 
band profile is considerably different. Empirical tight- 
binding calculations by Tanaka et al. gave a similar 
band structure along the F-M line, while the F-H dis- 
persion was not fully represented by the simplified model 
Hamiltonian. In a different paper, Tanaka et al. also 
pointed out that the dispersion of the conduction band 
is very sensitive to the oricntational configuration of the 
chains This is probably why the discrepancy with 
the DFT-LDA calculation for RbCeo by Erwin et al. |] 
is much worse: they used a model structure in which all 
the chains have the same orientation so that the unit cell 
includes only one AiCgo unit, while in the Reitveld model 
101 and therefore in the resultant structure of the present 
calculation, the adjacent chains are rotated by 90° from 
each other. 

Next, we would argue the character of the conduction 
band. The lower six conduction bands seem to be cat- 
egorized from their profiles along F-M line as three sets 
of two pairs as follows (FIG. 1.): looking at the energy 
levels at the F-point, the lowest and the third lowest, 
the second and the fourth lowest, the fifth and the sixth 
lowest, can be recognized as pairs. 

Direct inspection of these wavefunction ensures the 
above categorization. We have examined the iso-surfaces 
of these six wavefunctions at the F-point. Those belong- 
ing to the lowest and the third lowest are displayed in 
FIG. 2 (a) and (b), respectively. 

We can clearly see the bonding and the anti-bonding 
characters of the wavefunctions between the polymer 
chains in FIG. 2. (a) and (b), respectively. Furthermore, 
we can see the tiu character surviving in these states, i.e. 
the n* character of the C=C double bonds and the nodal 
plane which includes the four-membered ring connecting 
the fuUerenes. There is almost no evidence of hybridiza- 
tion with the other states. The ti^ symmetry and the 
bonding/anti-bonding characters as described above are 
also seen in the other four eigenstates, though the tiu 
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mark in the sixth eigenstate is weakened; the second (see 
FIG. 3.) and the fourth lowest eigenstates at the F-point 
have a nodal plane normal to the chain, while the fifth 
(see FIG. 4.) and the sixth lowest eigenstates have a 
nodal plane normal to the above two nodal planes. 

We can now understand that the wavefunctions of the 
conduction bands of the ortho-KCeo polymer are, at least 
at the F-point, linear combinations of the tiu molecular 
orbitals (MOs) as follows, 

*P,,r = ■^('/>p, ± 0pj, i = x,y,z. (0.1) 

Here, 4>p. and (pp. denote the MOs belonging to the first 
and the second fuUerenes in the unit cell, respectively. 
Note that the x-axis for the first fuUerene is taken to be 
parallel to the y-axis for the second fuUerene, since the 
coordinate axes of these p-like states for each fuUerene 
have been taken independently so that the x-axis includes 
the four-membered ring connecting the fuUerenes. 

Due to the polymerization, the threefold degenerate 
tiu states split into three different levels, Ep_^, Ep^, 
and Ep^ . These levels then constitute the inter-chain 
bonding/anti-bonding states, and the lowest state which 
has the inter-chain bonding character will be occupied. 
Clear inter-chain bonding character can be seen in FIG. 
2. (a), where the inter-chain interaction is large enough 
to form a conduction band with three-dimensional dis- 
persion. 

We note here that the occupied levels associated with 
the four-membered ring connecting the fuUerenes sink 
below the energy gap and so have not appeared in the 
above discussions. This may be confounding because al- 
kali doping is known to be crucial to the synthesis of the 
ortho-ACgQ. Below, we argue as to how the formation of 
the covalent connections couples with the profile of the 
conduction band. 

We now focus on the change in the local elec- 
tronic states around the carbon atoms forming the four- 
membered ring. Due to the polymerization, the local 
electronic states are converted from sp^ to sp^. This 
means that the a* and n* states of the sp^ re-hybridize 
into the a* states of the sp^. We emphasize here that (j)p^ 
has an amplitude around the ring (see FIG. 4), while (pp 
and (j)p^ do not. This means that only the levels which are 
functions of (pp^ will be affected by the re-hybridization. 
As can be seen from FIG. 1., the Px band is hfted about 
0.7eV from the Py-pz band, and forms a continuous band 
with the higher energy states. 

To confirm above the scenario in which the re- 
hybridization lifts only the Ep^ level, we have performed 
first-principles MO calculations on two different Cgo 
molecules having geometry equivalent to the polymer us- 
ing Gaussian 94 ||l^]. One is the bare Ceo molecule and 
the other is the Cgo molecule with the two carbon atoms 
which belong to the four-membered ring terminated with 
hydrogen bonds. In these calculations, Becke's 3 parame- 
ter hybrid method Jll| with Perdew's gradient-corrected 
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correlation functional jlj] (B3P86) and a 6-31G basis set 
were used. The bare deformed fuUerene shows a too- 
low Ep^ level, with Ep^ = -4.6eV < Ep^ = -4.2eV < 
Ep^ — —4. ley, while hydrogen termination serves to lift 
the Ep^ level, giving levels comparable to the band cal- 
culation, with Ep^ = -4.0eV < Ep^ = -SJeF < Ep^ = 
—2.9eV. Thus, we conclude that the formation of the 
four-membered ring connecting the fuUerenes raises the 
band from the continuous Py-pz band. In the latter 
MO calculation, the sequence, Ep^ < Ep^, is opposite 
from the F-point levels of our band calculation, but the 
difference is small in both cases. This is due to the crys- 
tal field effect, which is obserbed as the difference in the 
strength of the inter- fullerene interaction; in ^'p r, intra- 
chain as well as inter-chain bonding character can be seen 
between the fuUerenes, while the no intra-chain bonding 
character can be seen in 5'p^,r (see FIG. 3). 

Finally, we argue as to how the orientational configu- 
ration of the chains could be restricted as was reported 
experimentally There may be two key factors, one of 
which is the symmetry of the tiu molecular orbitals, the 
other being the rise of the Ep^ level above Ep^ and Ep^ 
as seen in our results. 

Nikolaev et al. calculated the tiu level splitting in the 
FCC AiCgo crystal as the function of the orienta- 
tional configuration of the fuUerenes by model calcula- 
tions in which the symmetry of the tiu state as well as the 
charge state of the fullerene are taken into account. They 
found that the two configurations which correspond to 
the local maxima in the energy splittings can be extended 
to represent the ortho-AiCeo (A=,K,Rb) and monoclinic 
A2A'C6o (A=Na,A'=K,Rb) polymers. Tanaka et al., us- 
ing the empirical tight binding method found that 
even in the orthorhombic AiCgo polymer, two orienta- 
tional configurations of the chains give local maxima in 
the band width. 

Although the threefold degenerate tiu eigenf unctions 
have p-like symmetries, they cannot be represented by 
spherical harmonics with £ — \. Owing to the icosahedral 
symmetry of the fullerene, they must be represented as 
the linear combinations of Y^m with £ = 5 ||lj| . The above 
studies suggest that the symmetry of the ti„ state is the 
major factor in determining the geometry of the ortho- 
AiCgo polymer. We then have to question as to why, if 
the model calculations show comparable energies for both 
local minima, AiCeo chooses an angle of 45 degrees. 

We can speculate that the asymmetry between the x- 
direction and the y-direction around the chain causes 
the energy difference between the orientational config- 
urations. The considerable rise of the px band from the 
continuous Py-pz band, or the difference in the diameters 
along the x-axis and along the y-axis may be the reason 
for the preference of the observed orientational configu- 
ration 

Further first-principles calculations of other orienta- 
tional configurations, or of polymers with different charge 
states, are required to clarify the role of the conduction 
electrons in determining the polymer geometry. 
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In conclusion, we have performed semi-SCF DFT-LDA 
band calculations as well as first-principles geometry op- 
timization of the ortho-KCgo polymer. The simulations 
were performed using cell dimensions equivalent with ex- 
periment. The resulting band has a three-dimensional 
dispersion which has a considerably different profile from 
those previously calculated ||^. We found that the 
lower six conduction eigenstates at the F-point are repre- 
sented as the linear combination of the tiu MOs. The rise 
of the Px band is found to be due to the local conversion 
of the electronic state, from sp^ to sp^, by the forma- 
tion of the four-membered ring connecting the fuUerenes. 
The conduction electrons will assist to form the one- 
dimensional chain as well as to bind the chains in the 
orthorhombic cell with this peculiar orientational config- 
uration. 
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Experiment This work 
Lattice parameter [A] 9.109(5) - 

9.953(5) 

14.321(5) 

Cl-Cl inter-fuUerene 1.65± 0.20 1.62 

Cl-Cl intra-fuUerene 1.74± 0.20 1.59 

C1-C2, C1-C3 1.50± 0.20 1.51 

TABLE L The calculated structural parameters and the 
experimental values. The notation is the same as Stephens et 
al. § 
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FIG. 1. Band structure of the ortho-KCeo polymer. The 
F-H direction is parallel to the chain and the F-M direction 
is normal to the chain. The vortical axis denotes the energy 
in eV. The origin placed at the fermi level. Note: the fermi 
level is not determined self-consistently; it is evaluated as the 
energy where the filling becomes half of the full occupation 
number during preparation of the effective potential by the 
CG method. 

FIG. 2. The iso-surfaces of the lowest, (a), and the third 
lowest, (b), eigenstate wavefunctions at the F-point are drawn 
with the atoms represented by a ball and stick model in the 
orthorhombic unit cell. The yellow balls denote the carbon 
atoms, and the white balls denote the potassium atoms. The 
viewpoint is parallel to the chain direction. The iso-level is set 
to 1/5 of the maximum absolute amplitude. The yellow and 
the blue curved surfaces denote the iso-surfaccs of the positive 
and negative values of the wavefunction, respectively. 

FIG. 3. The iso-surfaces of the second lowest eigenstate 
wavefunction at the F-point is drawn with the atoms repre- 
sented by a ball and stick model in the orthorhombic unit cell. 
The viewpoint is normal to the chain direction. 

FIG. 4. The contour plot of the fifth lowest eigenstate 
wavefunction at the F-point is drawn with the atoms rep- 
resented by a ball and stick model in the orthorhombic unit 
cell. The blue and red areas correspond to the negative and 
positive values of the wavefunction, respectively. The view- 
point is normal to the chain direction. The contour plane 
includes the four-membered ring connecting the fullerenes. 
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